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C D ■ Although consensus seems to exist about the validity of equations accounting for 

radiation reaction in curved space-time, their previous derivations were criticized 
C D ' recently as not fully satisfactory: some ambiguities were noticed in the procedure 

, of integration of the field momentum over the tube surrounding the world-line. To 

C ~ D ■ avoid these problems we suggest a purely local derivation dealing with the field 

quantities defined only on the world-line. We consider point particle interacting 
with scalar, vector (electromagnetic) and linearized gravitational fields in the (gen- 
erally non-vacuum) curved space-time. To properly renormalize the self-action in 
?H i the gravitational case, we use a manifestly reparameterization-invariant formula- 

tion of the theory. Scalar and vector divergences are shown to cancel for a certain 
ratio of the corresponding charges. We also report on a modest progress in extend- 
ing the results for the gravitational radiation reaction to the case of non-vacuum 
■ background. 

.=3. 

1. Introduction 

Study of the radiation reaction problem in classical electrodynamics initi- 
ated by Lorentz and Abraham by the end of the 19-th century , remained 
an area of active research during the whole 20-th century. Although with 
the development of quantum electrodynamics this problem became some- 
what academic, it still attracts attention in connection with new applica- 
tions and new ideas in fundamental theory. Current understanding of the 
radiation reaction has emerged in the classical works of Dirac 2 , Ivanenko 
and Sokolov 3 , Rohrlich 4 , Teitelboim 5 and some other (for a recent dis- 
cussion see 6 ). 

The Lorentz-Dirac equation was covariantly generalized to curved space- 
time with arbitrary metric by DeWitt and Brchmc in 1960. In their paper 
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( 7 ) an elegant technique of covariant expansion of two-point tensor quanti- 
ties was introduced which later became a basis of the perturbative quantum 
field theory calculations in curved space-time. It was used for calculation 
of the field momentum within a small tube surrounding the world-line of 
a point charge which resulted in the charge equation of motion with the 
radiation damping term. The main difference with the flat space case was 
the presence of the tail term depending on the entire history of a charge. 
Its presence signals violation of the Huygens' principle in curved space: a 
sharp pulse of light does not in general remain sharp, but gradually devel- 
ops a "tail". The equation (as extended by Hobbs to non- vacuum metrics) 
8 ) reads 

mz a = eF ia pf* + |e 2 (z Q - z 2 z a ) 

+ j^R^ + RyffP^z'^j +e 2 z p J' r 01 P(T')dr'. (1) 

where F m ^ is the incoming electromagnetic field R a p is the Ricci tensor 
and f a f3l is some two-point function taken on the world-line. This results 
was later generalized to other massless fields — scalar and linearized grav- 
itational. The gravitational radiation reaction was discussed long ago by 
C. DeWitt-Morette and Jing ( 9 ), and more recently reconsidered in detail 
by Mino, Sasaki and Tanaka ( 10 ) for vacuum background metrics. In 10 
an extension of DcWitt-Brehmc method was used consisting in integration 
of the gravitational field momentum flux across the small world-tube sur- 
rounding the particle world-line. The equation of motion obtained in ( 10 ) 
can be presented as 

2" = -\W + WWhX ~ Kv)* HP - ( 2 ) 

where (contrary to the case of (1)) the motion is supposed to be geodesic, so 
the local higher-derivative terms vanishes within the linear approximation, 
and the reaction force is given entirely by the tail term 

h£i = AmG £ w (G^ i A-^M^^. A ^(r)^(r')i T (r')i'(r'). (3) 

Although there is consensus about the validity of these equations, their 
previous derivation is somewhat problematic. In fact, as was discussed 
recently 11 , all calculations involving integration of the field momentum lo- 
cated the small tube surrounding the particle world-line contain some yet 
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unsolved problems. One such problem consists in computing the contribu- 
tions of "caps" at the ends of the chosen tube segment: contributions of 
"caps" were rather conjectured than rigorously calculated. Another prob- 
lem constitutes the singular integral over the internal boundary of the tube 
which was simply discarded in these calculations with no clear justifica- 
tion. In addition, the usual mass-renormalization is not directly applicable 
in the gravitational case since, due to the equivalence principle, the mass 
does not enter at all into the geodesic equations (in 10 the mass parameter 
was actually reintroduced by hand). 

Several new derivations were suggested during past few years. One is 
based on the redefinition of the Green's functions of massless fields in curved 
space-time proposed by Detweiler and Whiting 12 . But this scheme involves 
an axiomatic assumption about the nature of the singular term and so it 
docs not help to solve the above problems. Another scheme was suggested 
by Quinn and Wald 13 under the name of an "axiomatic approach to ra- 
diation reaction" . This scheme makes use of some intuitive " comparison 
axioms" , which do not follow from the first principles either. A recent at- 
tempt by Sanchez and Poisson 11 to overcome the above difficulties in fact 
is tight to a particular model of an extended particle (a " dumbbell" model) . 
Thus, though the results obtained within several different approaches agree 
in the final form of non-divergent terms, their consistent derivation from 
the first principles is still lacking. Also, the derivation of the gravitation ra- 
diation reaction force for non- vacuum background metrics remains an open 
problem. 

Here we briefly report on the new derivation (more detailed version will 
be published elsewhere) of the radiation reaction in curved space-time which 
has an advantage to deal with the fields only on the world-line and not to 
involve the volume integrals over the world-tube at all. This provides a 
more economic calculation and at the same time removes objections raised 
in 11 . The problem of divergences is relocated to the definition of the 
delta function with the support lying on the boundary of the integration 
domain, but this is exactly the same problem which is encountered in the 
flat space case where the approbated prescription amounts to the point- 
splitting procedure. Thus the local method is free from ambiguities which 
were met in the previous curved space calculations and can be considered 
as an adequate solution of the radiation reaction problem in curved space- 
time from the first principles. In general, our final results are conformal 
with those derived previously; in addition we partly remove the restriction 
by vacuum metrics in the gravitational case. Our signature is ( !- + +)■ 
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2. Non-geodesic motion 

Here we consider the non-geodesic motion of point particle along the affincly 
parametrized world-line x^ — z m (t) interacting with the scalar 4> and vector 
A^ fields. The total action is the sum S = S p + Sf, where the world-line 
part reads 



Sp = -m / (1 + q<f>)V^PdT + e A^z^dr, (4) 



(too being the bare mass, q, e - the scalar and electric charges) while the 
volume part is 



St = ~J (^-d^d^+^F^jJ—gd^x. 



(5) 



Our local approach to radiation reaction simply consists in the substi- 
tution of the proper fields into the particle equation of motion (for brevity 
we omit the external force) 

to (1 + qflzi* = eF^z" - m qU^ >. „ (6) 

where the velocity-transverse projector in the gauge i 2 = — 1 reads 

IF" = 0"" + 

The equations for the scalar field and the 4-potential (F = dA) read 

a<f> = 4np (7) 
aA^-R^A v = -4nj^ (8) 

where the covariant D'Alembert operator is understood □ = D^D^, and 
the source terms are standard. The Green's functions are defined in the 
usual way starting with the Hadamard solution (in the scalar case) : 



+ v In a + w 



(9) 



where a(x, z) is the Synge's world function, A is van Vleck determinant, 
and v, w satisfy the system 

Uv = 0, 

2v + (2v ;ft - ^)A- 1 A ;Jl^ )c7 ;^, + DA 1/2 + aUw = 0. 
In the vector case 

G { % = + v, a In a + w, a ) , (10) 
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where u fJ-a (x, z), v fJ-a (x, z) and w^a(x, z) are bi- vectors, satisfying a similar 
system. In particular, u^ a — j^A 1 ' 2 , where g„ a is the bi-vector of paral- 
lel transport. Notation is that of DeWitt and Brehme: indices a, (3, ... are 
associated with the "emission" point z, while fi, v, ... — with the "observa- 
tion" point x. When both points are taken on the world-line, we use the 
first set to denote z a (r') associated with the integration variable t', and 
the second set to denote z^(t), where r is the proper time in the resulting 
equation. In terms of these quantities the retarded solutions of the field 
equation read 

Trot (x) 

rct (x) =m q [ [-uS(a) + v6{-a)] dr' , (11) 



TrctO) 

F ™Ji x ) = ~ 2e J [u[^u]S' (a) + (u[ MQ;l ,] + v [fJia a u] )5(a) +v [va .^ z a dr' , 

— OO 

(12) 

where an anti-symmetrization over the indices /i and v is used. 

When substituted into the equations of motion, two-point functions 
become localized on the world-line, so we are led to use the expansions of 
the Synge's function <j(z(t), z(r'))(in terms containing delta-function and 
its derivative) as follows 

a(z(r), z(r')) = £ -— a(r,r)(r - r' f (13) 

k=0 

In what follows we will denote by dots the quantities 

a = (T a z a , a = a a fiZ a z^ + cr a z a , etc. (14) 

Denoting the difference as s = r — r', we get for a an expansion similar to 
that in the flat space 

^)--y-^ 2 M^ + 0( S 5 ), (15) 

but with dots corresponding to covariant derivatives along the world-line. 
Similarly, for the derivative of a with respect to z M (r) one finds 

S . ...„S Z \ . 4 



a»(s) = s ^i" - z*- + z* T + 0{s 4 ). (16) 

This quantity is a vector at the point z(t) and a scalar at the point z(t') 
where the index [i now corresponds to the point z(t): <t^ — da(z, z')/dz^. 
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An expansion for the delta-function reads 



6{-a) = S(s 2 /2) + s 



24 



5'{s 2 /2) + . 



(17) 



and an analogous expansion is easily obtained for its derivative. Since the 
most singular term is u5'(a)a v the maximal expansion order to be retained 
is s 3 . This also means that we have to retain for the delta- function only the 
leading term: 6(a) = 6(s 2 /2) + 0(s 4 ). The delta- function of the squared 
argument can be regularized in exactly the same way as in the flat space as 

d(s)5(s 2 /2) -» 6(s)8([s 2 - e 2 }/2) = 5{s - s)/e, 

where the positive regularization parameter e — > of the dimension of 
length is introduced. One finds in particular 



1 



es 



5'{s - e) 



. We have to perform the expansions up to the third order in s in all terms 
containing delta- function and its derivative (local terms). For the function 
u we will have Ricci-terms involved into expansions, up to third order one 
finds: 



u= 1 + ^R aT z T z a s 2 , 

and similarly in the vector case 

r 2 



ii-y _FL VT z s, 





a^z v z a = -sz^ + yz^ - s 3 Qz" + ^i? A „i A i^ + ^-i 2 ^ 



(18) 



(19) 



The functions of w-type in the tail term can not be found in a closed form. 

Collecting all the contributions, one obtains the following expressions 
for the scalar and electromagnetic reaction forces: 



/£ = mV 



(20) 



P =e 2 

Jem ° 



Divergent terms can be absorbed by the renormalization of mass. 

1 



m = too + 2^(m q - e ) 



(21) 



(22) 
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Note different signs of the divergent terms for scalar and vector self-forces, 
so for a "BPS" particle with the ratio of charges m\q\ = \e\ the model is free 
of divergencies. Finite contributions coincide with the previous results ob- 
tained via (somewhat questionable) world-tube derivations. Thus, in view 
of the criticisms in , our derivation may be considered as confirmation 
of these results. Technically, the local calculation is substantially simpler 
than the world-tube one. 



3. Neutral particle at geodesic motion 

In the case of gravitational radiation reaction we have the only one pa- 
rameter — particle mass, which actually does not enter into the geodesic 
equation. Thus it is not possible to use the above renormalization scheme. 
This problem can be remedied using a manifestly reparametrization invari- 
ant treatment. This is done by introducing the einbein e(r) on the world 
line acting as a Lagrange multiplier: 

™2 



S[z*,e]=-\ J 



e(r) g^z^z v 



e(r) 



dr, 



Varying (23) with respect to z M (A) and e(r) gives 
D 



dr 



(ei^) = 0, e : 



(23) 



(24) 



and we obtain the geodesic equation in a manifestly reparametrization in- 
variant form 



= 0. 



(25) 



yf-g^z^z" t 

We split the total metric into background and perturbation due to point 
particle g^ — ► g^ u = g^ v + nh^ v . Assuming now that the particle motion 
with no account for radiation reaction is geodesic on the background metric, 
the perturbed equation in the leading order in k will read 



uv _ 



{hx P - v - 2h uX . p )z A z p , 



(26) 



where contractions are with the background metric. 

The derivation of the field equations for the metric perturbation in the 
general case of non- vacuum background is non-trivial 14 . It is expected that 

(P) 

the particle stress-tensor T pv has to be divergence-free with respect to the 
background metric, this allows one to derive the equation for the metric 
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perturbation in a harmonic gauge, which leads to a manageable Green's 
function. But expanding the Bianchi identities for the full metric one finds 
that this is only possible if the Einstein tensor G M „ and metric perturbation 
hfj, v satisfies an additional equation 

G*h.x-Gx p h x i , <» = 0(z»). (27) 

Otherwise, the "naive" particle stress tensor on a given background is not 
enough, and construction of a reliable source term for metric perturbation 
becomes a complicated problem. Clearly, there is no general solution to the 
constraint (27), but some particular cases can be found. One can notice 
that this equation is identically satisfied for Einstein metrics = — A<? M „. 
Another case is conformally-flat metrics with some special conformal factor 
(details will be given elsewhere). Provided the Eq. (27) holds, one can 
derive (using the results of Sciama et al. 14 ) the following equation for the 
trace-reversed perturbation tp^ = — 

(P) 

ip^''^ + 2R^ v p ^P - 2ip^R^ a + ip^R - g^R a f3^ af3 = -2k 2 T^. (28) 

The corresponding Hadamard's function is similar to (10) with four-index 
bi-tensors u\ pa fj, v\ pa p. The retarded solutions is constructed in an anal- 
ogous way and substituted into the equation of motion (26) leading to 
an integral involving bi-tensor quantities depending on two points on the 
world-line r, r'. As before, one can distinguish four different contributions: 
terms proportional to the delta function, the derivative of the delta func- 
tion, the derivative of the Heaviside function, and the tail term. Performing 
series expansions in the first three cases up to the third order in r — r' = s 
one finds a local contribution to the self-force. To facilitate expansions of 
bi-tensors one has first to reduce them to scalars at the point of expansion, 
e.g. quantities like u v \ a j3Z OL (r')z !i (r') behave as a scalar with respect to the 
point z(t') and tensors with respect to z(t). Then expansion in s around 
r is straightforward. 

Contrary to the non-geodesic case, now we have to drop all (covariant) 
derivatives with respect to the proper time of the second and higher order, 
with the only exception of divergent term. Actually, the divergent term 
is proportional to the acceleration, so formally it vanishes for the geodesic 
motion. But still it may be argued that its infinite value demands some 
renormalization to be made, and we have to find an appropriate quantity 
to be renormalizcd . Since the mass does not enter the equation, we in- 
voke the einbein e(r), and insert as the corresponding background quantity 
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e =const. Collecting the expansions we obtain 



tail 



(29) 



Recall, that Ricci term is not arbitrary here, but has to be proportional to 
the metric for consistency. So actually the only local term here is propor- 
tional to the four-velocity z v and vanishes by virtue of the projector. The 
tail term is 



T 

" J [4v v \ a 0; p -2(9i>\vZ a p. p + V X pal3;u)- 



-gx p v: afi:v )]z' a z'Pz x zPdT'. (30) 
Renormalization of the einbein is performed as 



eo 



2e 



Finally the choice of the affine parameter z 2 
e = m, leads to the final form 



(31) 

-1, equivalent to setting 



z = mK 2 W v 



[(2lVAa/3;p - gv\g aT V aTa ^p) 



- {vx paP :u - 1/2 gXpg^Vvrafr^jz^z'PdT'^zP] 



(32) 



This coincides with the result obtained in 10 > 13 . We have thus shown that 
this equation remains valid for a class on non- vacuum metrics, in particular, 
for Einstein spaces. 



4. Discussion 

We have presented a local calculation of radiation reaction in both geodesic 
and non-geodesic cases which do not involve any ambiguous integrals out- 
side the world-line and thus is free from the corresponding problems. It 
is based on the expansion of bi-tensor quantities only on the world line 
and technically is much simpler than the DeWitt-Brehme type approach. 
We have also generalized the gravitational radiation reaction to some non- 
vacuum metrics, however, satisfying an additional condition. The general 
non- vacuum case can not be treated within the test body/external field 
approximation because the global Bianchi identity demands to take into 
account the perturbation of the background matter as well. 
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